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We study the role of quantum statistics in the performance of Otto cycles. First, we show 
analytically that the work distributions for bosonic and fermionic working fluids are identical for 
cycles driven by harmonic trapping potentials. Subsequently, in the case of non-harmonic potentials, 
we find that the interplay between different energy level spacings and particle statistics strongly 
affects the performances of the engine cycle. To demonstrate this, we examine three trapping 
potentials which induce different (single particle) energy level spacings: monotonically decreasing 
with the level number, monotonically increasing, and the case in which the level spacing does not 
vary monotonically. 


I. INTRODUCTION 


A major technological challenge today is the design of 
devices at the nanoscale that can efficiently and reliably 
convert energy between different forms. At such small 
scales, the thermal and quantum fluctuations of thermo¬ 
dynamic quantities become especially significant, and the 
stochastic nature plays a dominant role in determining 
the performance of quantum devices [IHl- 

More recently, major theoretical developments have 
also been made specifically towards the understanding 
and implementation of quantum thermodynamics. These 
include for instance, advances in fluctuations theorems 
[3) H) [1| , and on the use of external controlprotocols to 
produce adiabatic dynamics in finite time |7Hl4l|. 

At the same time, properties of non-equilibrium work 
distributions have also been investigated experimentally 
in classical systems [T^ [l^ , culminating with the realiza¬ 
tion of a colloidal micro-scale stochastic heat engine (l3|. 
and quantum systems [l^. Due to the high degree of 
tunability and control, ultra-cold ion setups have also 
positioned themselves as promising candidates for exper- 
imental realizations of quantum heat engine cycles |19j |. 
In addition, a proposal for the realization of an Otto cycle 
within a solid state set-up was recently studied [^ . 

More recently, it has also been shown numerically, that 
the exchange symmetry of particles do play a significant 
role in the work distribution of quantum many-body sys¬ 
tems [llj . As such, here we address the pertinent issue of 
many-body statistics in quantum heat engines. In par¬ 
ticular, we elucidate the role of exchange statistics in the 
performance of quantum cycles by providing an intuitive 
understanding of the interesting connection between ex¬ 
change statistics and the geometry of the trapping po¬ 
tential. 

We consider a quantum equivalent of the Otto-cycle 
in which the working fluid is a quantum gas of N 
indistinguishable and non-interacting particles. The col¬ 
lective system undergoes a cyclic sequence of adiabatic 
and iso-parametric processes, where the heat received by 
external reservoirs can be in part, converted into useful 
work, while the rest is transferred to a cold reservoir. Im¬ 
portant characteristics and figures of merit of the engine 
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FIG. 1. (Color online) (a-c) Different potential wells with 
corresponding single particle energy spectrum Sj. The dark 
blue (light red) squares represent the occupation of energy 
levels at 0 temperature for the case of 4 bosonic (fermionic) 
atoms, (a) triangular potential (o = 1), (b) infinite square 
well potential (a = oo) and (c) double well potential. The 
horizontal thin grey lines indicate the position of the j-th 
single particle energy level, (d) Schematics of a quantum 
Otto cycle in an average energy (E) vs trapping parameter 

OJa- 


cycle’s performances are the probability distributions of 
both its net work output and its efficiency. In this work 
we study how the statistics, either bosonic or fermionic, 
of the gas’ particles affects the engine’s performance. We 
show that this is strongly dependent on the single particle 
energy level structure of the system and indicate strate¬ 
gies to improve the performance of the engine cycle by a 
careful choice of both the trapping potential and the type 
of gas. It is only when considered in conjunction, that 
trap geometry [2^ and particle exchange statistics [2l| 
can be utilized to optimize the performance of quantum 
heat engines. 
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We consider non-interacting particles in different trap¬ 
ping potentials: with the energy level spacings either 
constant (harmonic potential), monotonically decreasing 
(triangular potential, see Figllja)), increasing (square 
potential, see FigITJb)) or change non-monotonically 
(double well, see FiglTJc)). To gain some preliminary 
intuitive insights to the problem, in FiglT] we show how 
4 atoms are distributed in the respective potentials at 
temperature T = 0 depending on the given symmetry. 
We note that the energy difference between the ground 
state and the first excited state can be very different de¬ 
pending on whether the gas is bosonic or fermionic in 
nature and on whether the gap between energy levels is 
increasing (e.g. square well) or decreasing (e.g. triangu¬ 
lar potential). In FigUJa) we observe that for bosons 
(blue squares) the gap is larger than for that of the 
fermions (red circles); Whereas in the case of a poten¬ 
tial in which the energy level distance is monotonically 
increasing (square well in FiglTJb), the gap will be larger 
for the fermions and smaller for that of the bosons. For 
a potential in which the distances between energy lev¬ 
els do not feature a monotonous behavior (double well 
in FigHKc)), the excitation gap not only depends on the 
particle statistics, but also on the number of atoms. 

This paper is organized as follows: In section |IT1 we 
provide general relations and results for quantum Otto 
cycles with scaling potentials; In section nil we show 
our main results and lastly, in section m we draw our 
conclusions. 


II. THE OTTO CYCLE IN SCALING 
POTENTIALS 

In classical thermodynamics, the Otto cycle is com¬ 
posed of two processes in which only heat is exchanged 
with the environment, and two processes in which only 
work is exchanged (and no heat is transferred). The for¬ 
mer is usually referred to as an isochoric process for a 
classical ideal gas. Analogously, a quantum Otto cycle 
is composed of two quantum adiabatic processes, repre¬ 
sented in FiglTJd) by the steps [A] [B\ and [C] [D\, 

as well as two iso-parametric processes, represented in 
FigUDd) by the steps \B] —>• [C] and [D] —>• [A], in which 
the Hamiltonian of the system is kept constant and the 
environment is (weakly) coupled to the system . Sim¬ 
ilar to their classical counterparts, no heat is exchanged 
with the working fluid in quantum adiabatic processes, 
while only heat is exchanged in an iso-parametric pro¬ 
cess. T[^] corresponds to the lowest temperature in the 
cycle while T[c] is the largest. As we are interested in 
the governing principles of quantum thermodynamic cy¬ 
cles involving non-interacting many-body quantum gases, 
we will only study quasi-static cycles, leaving the cycles 
with finite power for a future study. 

To gain deeper analytical insights, we also focus on 
Hamiltonians with scaling properties [H, [H, [2^ , or 

more generally on driving protocols in which the energy 


level structure after a unitary process is simply multi¬ 
plied by a scale-factor. The harmonic oscillator is the 
prototypical example of such scaling Hamiltonians, as the 
energy levels are simply proportional to the trapping fre¬ 
quency. 

In a more general sense, for every potential of the form 
(where A is a control parameter with 
units of length and Xi is the position coordinate of the 
f-th atom), the Hamiltonian H has a scaling energy spec¬ 
trum. In fact 
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where = x^/A is a dimensionless variable [l^. This 
applies to all power-law potentials of the form V{xi) = 
^m{uJa\xi\)°' since 
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where Ca(wa) = {hjxiaf°‘ l’m°‘ ^ is the energy scale 
[ 2 ^ [ 2 ^. It is thus now clear that that the energy levels 
are simply scaled by the prefactor Co(wa)- For instance, 
the ratio of the energy of the j-th level for Hamiltonian 
parameter a;(, (i.e. ej(a;(,)) divided by that for parameter 
Wa (i-C. ej{u!a)) does not depend on j and it is given by 

Sji^a) (a(x’a) 


Ca(^o) 


( 3 ) 


where ^ is a parameter that only depends on Wa and 
In the following we will use the lighter notation Ca(<j-'a) = 
Ca and Ca(Wa) = Ca¬ 
in addition, a more general polynomial like the double¬ 
well potential may also be made to possess a scaling prop¬ 
erty in its energy levels, albeit in the presence of certain 
restrictions. For example, the double well potential 

y{xi) = '^ {-ujIx] + ujjxt) (4) 


has a scaling Hamiltonian 


h wl 4 
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( 5 ) 


so long as the anharmonic parameter 7 = huif/mui^ re¬ 
mains fixed throughout the Hamiltonian evolution. 


A. Mean Work {Wp) 

Let us consider a process from Hamiltonian parameter 
LOa to Lo'^. The respective single particle j-th energy eigen¬ 
values will be Ej and e'. Let rij represent the number of 
particles in the j-th single energy eigenstate. 
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A single A^—body state is thus given by the ordered 
set of occupation numbers {rij}. A process is performed 
quasi-statically by changing the trap parameter adiabat- 
ically, and the corresponding TV—body microstate at the 
end of the compression is then given by {m-fc}. For a sin¬ 
gle realization of the protocol p the work is given by a 
two-times energy measurement [3, Hi 

Wp = - Elm,} ( 6 ) 

The energy (resp. given by the sum 

over the energy si (resp. e() of the Z-th single particle 
eigenstate of the Hamiltonian (e.g. = J^k^knT-k)- 

Hence, the ensemble average of work done for the process 
between two trap parameters is given by: 

(Wp) = "YL ^ ^E{m,} ~ E{nj})^{m,},{nj}P{nj} 

{rij} {m,} 

= Y.{eU^}-Em)Pm (7) 

{"jl 


where the Kronecker delta S{m,},{nj} is due to the adi- 
abaticity of the process. The average (...) is computed 
over the intial probability distribution P{nj} which, at 
this stage, can be general. Because of the scaling prop¬ 
erty of the Hamiltonians we consider, the ratio of ini¬ 
tial to final energy of each single particle eigenstate 
{Sjjsj) = V (see Eq.(l3]) and [H, H^). This results in 
the work of an adiabatic process to be 


It follows that the scaled standard deviation of work, 
defined as the standard deviation of the work distribution 
divided by its mean work is given by 


_ / {w^) . 

{W,) - V {W,)^ V {EY 


( 11 ) 


which is independent of i.e. independent of whether (i) 
the process is a compression or an expansion, and (ii) the 
extent of the process (provided it is quantum adiabatic). 


C. Net work output of a cycle: average and 
variance 

The ensemble moments of the cyclic work distribution 
is derived from the work distribution of the constitutive 
single adibatic processes. The cycle work distribution is 
composed of two statistically independent work distri¬ 
butions on the compression and expansion arms of the 
cycle because of the nature of the ideal quantum Otto 
cycle. The combined work distribution is thus given by 
the product distribution of two statistically independent 
processes 

P{W) = (5[IT - {W, + We)]P{W,)P{We) (12) 

where Wc {We) is the work in the compression (expan¬ 
sion) process. For the cycle depicted in FiglTJd) between 
parameters Wa and the net work output is given by 


{Wp) = {y-l) {E) ( 8 ) 

where 

(E) = ^WjP{ni} (9) 

{"j} i 

is the energy of the system at the beginning of the pro¬ 
cess. 


B. Work Standard Deviation awp 

The standard deviation of the work probability distri¬ 
bution function of a process p, uwp = \J{W^) — {Wp)'^, 
is also important to characterize a process as it quanti¬ 
fies the consistency and reliability of the output. For an 
adiabatic unitary process in a scaling potential we have 

{Wp) = ^ ^ ~ E{nj}) ^{mk}{nj}P{nj} 

{rij} {m,} 

= (^- 1 )" I E{n,} 

W} V T / 

= {i^-1)Ue^) ( 10 ) 


{W)e, = {v- 1 ){E^A])c. + (^ - l) (13) 

where here again, v = Ca/Ca and (iT[y])a is the average 
energy at the vertex [V] of the cycle and (... )„ means 
average over the thermal canonical state for fermions 
{a = F) or bosons {a = B). 

Consequentially, the mean and variance of the cycle 
are also just the linear sum of the mean and variance of 
the individual processes. The scaled work fluctuations 
are thus given by 

^ ^/(AVF,J2 + (AWeJ2 

(bF)a {Wc)c. + {We)c. ^ ’ 

In addition, we also note that if the adiabatic processes of 
an Otto cycle are peformed within scaling potentials then 
the efficiency is a non-stochastic function of v inde¬ 
pendent of the quantum statistics of the working fluid. 
At each realization of the cycle 

W 

_ (1 ~ l/^)^[C],{mfc} + (1 - ^)E[A],{nP 
ElC],{m,} - l'ElA],{u^} 

1 

= 1 - 

V 


( 15 ) 
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where indeed v > 1 and Q[b]^[C] is the heat transferred 
into the system between states [B] and [C]. In Eg. dTSl) 
E[v],{mk} represents the energy of the eigenstate {mk\ 
at vertex [V] of the cycle. 


III. CYCLES PERFORMANCE FOR BOSONIC 
AND FERMIONIC GASES 

The average work output and its standard deviation 
can thus be computed from Eq.dH) and (ITO)) . The av¬ 
erage energy and variance of a canonical thermal state 
is computed using {E)a = —d{\nZ'^) /d(3 and = 
{In Zfi) /(3 where (3 = l/ksT is the inverse tem¬ 
perature and ks is the Boltzmann constant. To obtain 
accurate results for low number of particles and low tem¬ 
peratures, we use the canonical par ition function Z'^ for 
N (fermionic or bosonic) atoms 2^. It can be computed 
using the following recursion relation for non-degenerate 
spectra 


Z%{I3) = ^ ^(-l)^Z“(n/3)Z“_„(/3) (16) 

n—1 

where 7 = 2n for bosons (a = B) [ 2 ^, [ 2 ^ and 7 = 
n + 1 for fermions (a = F) [s^ and where ^“(/3) = 
is the single particle partition function (which 
is obviously independent of the value of a). Note also 
that Zq =1. Ea. (IT6)l (derived in detail in the Appendix 
A), will be used to compute Z'^{l3) numerically. Eor the 
harmonic case, exact analytical results will be described 
in the following section (III-A). 


A. Harmonic potential 


Eor a simple harmonic oscillator of frequency a; 2 , it is 
possible to compute the partition function with the use 
of the q—shifted factorials {z-,q)n = n^=i(l “ zq‘^~^). 
The partition function for N bosons is given by 


^-N0huJ2/2 

(<?; q)N 


and for N fermions by 


^-N^Phuj2/2 

(<?; q)N 


(17) 


(18) 


where q = (see Appendix B for a detailed deriva¬ 

tion). This results in average energies, (£!)„, given by 


{E)b 


fuV2N 

2 


N 




huj2n 

^0huj2n _ ^ 


(19) 


{E)f = {E)b + ^N{N-1). (20) 






Sa/kB Sa/kB 

FIG. 2. (Color online) Temperature T vs entropy Sc diagram 
for an Otto cycle in an infinite square well potential (a = 
00 ) with fermion (red-dashed line) or bosons (blue-continuous 
line) as working substances. The trapping parameter varies 
between ^00 and = 2(oo while the minimum temperature 
T^a] = 2 000 . The particle number N and the maximum 
temperature Tyc] takes the following values: (a) A = 2 and 
T[c] = 5 0OO, (b) A = 2 and Tyc] = 100 0oo, (c) A = 5 and 
Tycy = 5 0OC, (d) A = 5 and Tycy = 100 0oc. 



FIG. 3. (Color online) Temperature T vs entropy Sa dia¬ 
gram for an Otto cycle in a triangular potential (a = 1) with 
fermions (red-dashed line) or bosons (blue-continuous line) as 
working substances. The trapping parameter varies between 
Cl and Cl = 2Ci while the minimum temperature TyAy = 0i/3. 
The particle number A and the maximum temperature Tycy 
take the following values: (a) A = 2 and Tycy = 0i, (b) A = 2 
and Tycy = 10 0i, (c) A = 5 and Tycy = 0i, (d) A = 5 and 
Tycy = 10 01 . 


and 
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All higher moments of the energy are identical for bosons 
and fermions. In particular, the standard deviation for 
bosons geb is equal to that of fermions gep 


<^Eb = '^Ep 


N 


E 




( 21 ) 


The constant mean energy difference between fermions 
and bosons ensures that the average work transfer is al¬ 
ways greater for that of fermions, (bbp)_F > while 

the rescaled standard deviation is instead, always larger 
for bosons < ^w’’fB 0" course assumes 

that the comparisons are between given initial canonical 
states of the same temperatures and identical quantum 
adiabatic protocols in a harmonic potential. 

For the work output of the full cycle, with Hamiltonian 
parameter that changes between W 2 and W 21 using Eq. CSD 
one can also easily show that 


{W)f 


{W)b + 

{W)b + 
{W)b 


{v - 1)W2 -t- - 

' V 



A 



W2 -I- 




A 

( 22 ) 


where A = |A(A- 1). 

Hence for a reversible quantum Otto cycle driven by 
harmonic potentials with non-interacting particles, all 
moments of work distribution for the cycle are the same, 
regardless of the working fluid being bosonic or fermionic. 


B. Triangular and square well potentials 

The previous results for the harmonic oscillator are a 
direct consequence of the constant energy level spacings. 
For non-harmonic potentials, where the level spacing is 
not constant, the results will be very different. For in¬ 
stance, from Figllja) it is clear that, at zero temperature, 
a fermionic system requires less energy to be excited com¬ 
pared to a bosonic one (because the gap between consec¬ 
utive energy levels monotonically decreases, hence the 
energy gap from the last filled state to the first unfilled 
one is larger for bosons than for fermions). The converse 
is true for a square well potential in which the energy 
level spacing monotonically increases, see FiglTJb). As 
we will show later, this has strong consequences on the 
performance of an engine cycle. 

To describe this in a clear and pictorial manner, see 
Figsi2]and|31 we use a representation of a cycle typical of 
classical thermodynamics, the temperature vs entropy T- 
S diagram (where T is divided by a scale of temperature 
given by 0a = Ca/^B while in the cycle the Hamiltonian 
parameter varies between uja and > Wq). To compute 
the entropy Sa we use the thermodynamic relation = 
{E)a — TSa, where the free energy = —fcBTln(Z^). 
In this diagram, the area within the cycle represent both 
the net heat and net work exchanges. 



FIG. 4. (Color online) Fermionic to bosonic ratio of the mean 
work output of an Otto cycle vs. high temperature T^c] for 
different particle numbers in the following potentials: (a) Tri¬ 
angular potential (a = 1), (b) Infinite square well potential 
(a = 00 ) and (c) Double well potential (a = 2, 7 = 0.01). In 
all potentials, C,'^ = 2^. In (a) = 0i/3, while in (b) and 

(c) T^a] = 20a. Different numbers of atoms are represented 
by the different (colored) lines: N = 2 for the dashed green 
line, A = 3 for the dot-dashed red line, A = 4 for the dot¬ 
ted black line and A = 5 for the continuous blue line and 
A = 10 for the dotted brown line in (c). The insets depict 
the potentials used. 


Using Eq. (HU) we are able to compute Z'^ for different 
particle numbers and cycle parameters. 

We consider only cases in which the coldest temper¬ 
ature T[a] is cold enough to see the combined effect of 
both energy level spacings and particle statistics [s^. We 
then consider different cases in which the value of the 
highest temperature T^c] j varies from (relatively) cold to 
hot, and also with working fluids that consist of different 
number of particles. Eor low T^c] we expect large differ¬ 
ences in the distribution for bosonic or fermionic gases. 
In EigH] we depict the case of an infinite square well, 
while in EigO we study the case of a triangular poten¬ 
tial. In both figures the fermionic cycle is represented 
with dashed red lines and the bosonic cycle by continu¬ 
ous blue lines. Whether a bosonic or a fermionic working 
fluid produces more average work output depends on the 
difference between the adiabatic work in the cold (com¬ 
pression from [A] —^ [H]) and in the hot arm (expansion 
from [C] ^[D]). 

As clearly shown in EigHKc), when the number of par¬ 
ticles are large and both temperatures T^a] ^-nd T^c] 
cold enough, there can be remarkable differences between 
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FIG. 5. (Color online) Fermionic to bosonic ratio of the stan¬ 
dard deviation in the work distribution of an Otto cycle vs. 
high temperatures T[c\ for different particle numbers. Values 
for the cold temperature Tja] , particle number N and Hamil¬ 
tonian parameter values uJa are as given in Fig[4l The insets 
depict the potentials used. 


the work output of a fermionic or a bosonic system. The 
area enclosed in the fermionic cycle is much smaller. In 
fact, for fermions the energy level separation is too large 
for entropy to change significantly. A similar scenario, 
although in the opposite direction, emerges for a trian¬ 
gular potential as depicted in Figl3l)c). In this case the 
energy level spacing is larger for the bosons and the en¬ 
tropy will be lower (and less heat is absorbed). For small 
number of atoms and large enough highest temperatures, 
the ratio of the net work output (IF)f/(IF) b tends to 1, 
see FigsHKb) and|3Kb). 

The case of two atoms and cold temperatures are de¬ 
picted in FigsISKHKa). In this scenario the number of 
atoms is too small to result in a marked difference of the 
cycles due to the quantum statistics. Cycles for larger 
number of atoms and hotter temperatures are shown in 
FigsI^HSKd). It should be noted that for the triangular 
potential, FiglJKd), almost no difference can be noted 
between bosons and fermions in the high temperature 
regime. However, in the cold temperature portion of the 
cycle, important and large differences in the amount of 
entropy exchanged can be noticed. In particular, the en¬ 
tropy is lower for the species experiencing the largest gap 
(bosons in this case). 

To gain deeper understanding we revert to Figsl31[SKa- 
b). The triangular and the square well potential have 
qualitatively opposite behaviors. When the highest tem¬ 


perature Tjc] is small, the total entropy and the entropy 
change for the bosons in the triangular potential (re¬ 
spectively for the fermions in the square well potential) 
are much smaller than that of the fermions (respectively 
bosons), see Figsl^Kc) and [UJc). It is easier to under¬ 
stand the qualitative behavior for large T]^c] by consider¬ 
ing the work exchanges. In this case the ratio of the work 
exchanges {W[c]^[d])f/{W[c]^[d])b = {E[c]) f/{E[ c ])b 
and thus it will eventually converge towards 1 as the tem¬ 
perature increases. However the difference in the work 
processes {W\^c]^[d])f — {W\^c]^[d])b behaves differently 
depending on the type of potential. It increases indefi¬ 
nitely for the square well, while it goes to 0 for the trian¬ 
gular potential. Hence, for the square well potential the 
ratio {W)fI{W)b will tend to 1 from above, while for 
the triangular potential it will tend to 1 from below . 
The same analysis also explains the qualitatively similar 
behavior of the variance of the work probability distribu¬ 
tion, Figini For this quantity, however, the crossings of 
the line opjoB = 1 occur at a different values of T]^c]- 

C. Double well potentials 

The double well (Figlljc)), is a typical example of a 
system in which the energy level spacing does not vary 
monotonically. Hence, we will use this to study the in¬ 
terplay of geometry and particle statistics in the perfor¬ 
mance of the engine cycle. As shown earlier in section 
im the double well potential can have a scaling property 
so long as 7 = hujf/muj^ is a constant. For ease of ana¬ 
lytical and numerical analysis, we will focus on this case. 
As shown in FigHJc) the fermionic to bosonic ratio of 
the work output {Wf)/{Wb) has a pronounced oscilla¬ 
tory behavior compared to potentials with monotonically 
changing energy level structure as the (inverse) temper¬ 
ature of the hot reservoir f3h- 

The same behavior is expected and observed also for 
the ratio of the standard deviation of the work done 
for a fermionic over a bosonic working fluid CTi^’/o's, see 
FigEJc). 

IV. CONCLUSIONS 

In this work we have studied how quantum statistics 
of particles in a working fluid can affect the performance 
of an engine cycle. We have shown that the interplay 
of statistics with energy level structure of a system has 
remarkable and non-trivial effects on the engine cycle’s 
output. We have analyzed systems with monotonical- 
lly increasing (square well potential) and decreasing (tri¬ 
angular potential) energy level structure, which possess 
qualitatively opposite behaviors, both in the low and 
high temperature regimes. For potentials with a non¬ 
monotonic energy level spacing (for example, a double 
well potential), the engine cycle’s peformance is charac¬ 
terized by the presence of oscillations in the fermionic to 
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bosonic ratio of a given figure of merit, and the proper¬ 
ties of these oscillations are closely associated to the un¬ 
derlying energy level structure of the Hamiltonian, the 
particle number and the temperature regimes. As a final 
remark, systems driven in finite time, and many-body 
interactions will be the subject of future studies. 
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Appendix A: Derivation of recnrsive relation for the 
canonical partition function 

Here, we derive Eq. (1161) following the argument as pre¬ 
sented in [1^ . To do so we start from the grand canonical 
partition function where a = i? or E, is the subscript 
that denotes bosonic and fermionic statistics respectively. 
The grand partition function can be written as 

OO 

= 1 + ^ (Al) 

N=1 

where is the canonical partition function for N atoms. 
In equation m z = is commonly known as the 
fugacity, where ^ is the chemical potential associated to 
the ensemble. 

It is also possible to write the logarithm of Sq, as 

OO 

(A2) 

N=1 

where in general, assumes different functional forms 
depending on the given type of particle statistics con¬ 
sidered. By exponentiating (IA2p and using the Eaa di 
Bruno’s formula, which states that the power series ex¬ 
pansion of the exponential of a polynomial /N\ 

is 

exp I /N\ j = B{ai, ...aN)x^/N\ (A3) 

\ N / N=0 




we can derive 


Appendix B: Partition function for the simple 
harmonic oscillator 


= e(^”=i Bn (Ka,l> ■■■: ^ 

N^O 

(A4) 

where Bn is a Bell polynomial. Next, by equating corre¬ 
sponding terms of the polynomial expansion in (IA4I) and 
(TATI) , we get 


YCx _ 
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Bn^k, 


CK,! 5 


J 


N\ 


(AS) 


Note that up to this point in the derivation, we have 
not made any distinction between bosonic or fermionic 
statistics. The grand partition function for the fermions 
can be written as a formal power series 


ln(.=,ir) = ^ In (l-I-ze 

k 

k N 

= - ly. A 

N \ k / ' 

= J2{-l)^+\N-l)\Zy{Nf3)— (A6) 

N 


where we have used the expansion of ln(l -l-a;) and where 
Zi(n(3) is the single particle canonical partition function 
with inverse temperature n[3 (of course Zy = Z^ = Z^). 
Finally, using equations (IA5|) and (IA6I) together with the 
identity for Bell polynomials; 


For the simple harmonic oscillator, the Hamiltonian 
parameter a = 2. Hence the energy of the system will 
scale with huj 2 - Following the derivation of [^, we be¬ 
gin by writing the grand canonical partition function for 
bosons in terms of a q-shifted factorial: 


^Biz )= n — 


fe =0 


q'^z 


{z]q)c 


(Bl) 


Here, q = z = where z = is as 

before, the fugacity, and (z;g)oo is a q-shifted factorial. 
By using the g-analog binomial theorem: 


1 

(■5; q)co 


OO 


E 

n=0 



(B2) 


where the ^-shifted factorial reduces to simply: 


N 

x—1 


Notice that we can also write r2(z) as a formal power 
series involving the canonical partition function Z'^ 


H„(z) = 1 + ^ Z^z'^ (B4) 

n—0 

And thus by comparison of the two equivalent series, we 
are able to write down the canonical partition function 
for the bosons: 


(^1, A ^ ^ (771 1 ) 

m=l ^ ^ 

(A7) 

we derive Eg. (1161) for fermions. Similarly for bosons 
In(SB) = — y] In (1 — ze“^^*’) 

k 

____ 

k n 

= ^{n-l)\ Zi{nl3)— (A8) 


H, 


and again using Eo. dASI) . (IA7I) and (IA8I) yields Eo. lfTO 
for bosons. 
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( 9 ; q)N 

Similarly for the fermions, 

flpiz) = {-z;q)oo 
and using the identity 

^ ^ {-l)NqNiN-l)/2 ~^N 

{z,q)^ = 2^ — 


we can write 


N^O 
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(B6) 


(B7) 
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